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Numerical simulations of force-free, degenerate (ffde) pulsar and black hole 
magnetospheres are often based on 1-D characteristics. In particular, the plasma 
wave polarizations that can be propagated along the 1-D characteristics deter- 
>■ ■ mine the time evolution of the entire system. There are two sets of characteristics, 

corresponding to the fast and Alfven modes. The fast wave is generally consid- 
er ' ered to be a transverse light wave, however recently it has been claimed that 
Q ■ light-like fast waves can transport a longitudinal electric polarization, E\\., at the 
^ ■ speed of light. The implication is quite profound if true, namely that the wrong 
' information has been propagated along the fast characteristics in all previous 
Q^' simulations of force-free magnetospheres. It is shown in this Letter that the 
■ light-like fast waves must be transverse and previous simulations are valid. This 
result is demonstrated by means of a fundamental physical principle (associated 
c3 ■ with the fact that particles cannot flow faster than the speed of light), there exists 
a charge horizon in ffde magnetospheres. It is shown that the Alfven critical sur- 
^ ■ faces in a ffde magnetosphere are both charge and particle horizons, i.e. one way 
■ membranes that do not permit traversal by charges nor particles anti-directed to 
the bulk flow. Since the propagation of a discontinuous change in £"11 requires 
a physical surface charge on the wave-face, it is also a one-way membrane for 
longitudinally polarized waves. Besides justifying previous ffde simulations this 
result also invalidates previous claims that fast waves can radiate E\\ from the 
event horizon of a black hole. 



1. Introduction 

The zero mass, cold limit of MHD (magnetohydrodynamics) known as ffde is useful for 
studying tenuous plasmas in strong magnetic fields that might occur in pulsars, black holes 
or gamma ray bursts (Blandford 2002). Recently, great progress has been made in pulsar 
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ffde simulations, Spitkovsky (2003), and in ffde simulations of black hole magnetospheres, 
Komissarov (2001, 2004); Uzdensky (2003). To varying degrees each of these simulations 
require knowledge of the 1-D characteristics of the ffde system in order to time evolve the 
magnetosphere. Specifically, the polarization properties of the waves determine the changes 
in the fields that can be propagated at the appropriate speed along a particular characteristic 
direction. The 1-D characteristics represent local plane wave (no transverse variation of 
the field values along the wave face) solutions that can be pieced together on a mesh to 
produced higher dimensional wave structures. Thus, it is crucial that the 1-D characteristics 
and the corresponding local plane waves are well determined before utilizing the method of 
characteristics to time evolve a magnetosphere. There are two plasma modes in the system, 
the fast mode and the Alfven mode. It has been previously shown that a light-like fast wave 
has a transverse electric polarization (Blandford 2002; Uchida 1997; Komissarov 2002). 
However, recently in Levinson (2004), it was claimed that light-like plane wave solutions 
can attain a longitudinal polarization (k ■ E 7^ 0, where k is the propagation vector) in a 
rotating magnetosphere, a case previously considered in Komissarov (2002). The polarization 
information of Komissarov (2002) was used to piece together 2-D simulations of black hole 
magnetospheres by means of a Riemann solver (Komissarov 2001, 2004). If the the claim 
of Levinson (2004) is correct then the simulations of Komissarov (2001, 2004) are based 
on the wrong information being propagated along the characteristics. This Letter tries to 
resolve this controversy by studying the global physical properties of a ffde magnetosphere 
as opposed to debating over mathematical errors. In particular, we find that the Alfven 
critical surface is a one-way membrane for the fiow of charge, a charge horizon. A similar 
result was previously shown in Komissarov (2004). It is demonstrated that the propagation 
of a discontinuity in £"11 requires a physical surface charge on the wave-face by Gauss' Law. 
Hence discontinuities in £"11 (step waves) cannot traverse the charge horizon, anti-directed 
to the flow. It also precludes the possibility of information on propagating away from 
the event horizon. Finally, the association of with the Alfven wave yields a nice physical 
interpretation of the Alfven critical surface as the charge horizon. 

The Letter is structured as follows. First of all, the maximum velocity of a particle is 
determined in ffde in section 2. With this concept, it is straightforward to deflne the charge 
horizon in section 3, with an example in a black hole magnetosphere presented in section 4. 
In section 5, it is shown mathematically that any propagating discontinuity in E\\ must travel 
at the Alfven wave speed, thus it cannot travel isotropically at the speed of light. Physically, 
by Gauss' law the wavefront is a propagating surface charge, hence it cannot cross the charge 
horizon anti-directed to the flow. 
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2. Force- Free Magnetospheric Dynamics 

The forcc-frcc conditions are given by the following relationships, written covariantly in 
terms of the Maxwell field strength tensor, F^^ , and four-current density, J^, as well as in 
component form, 

F'^'^J.^O, p,E+^^ = 0, (2-1) 

c 

with the implied constraints, 

F^^^F^, > , *F'^'^F^, - . (2-2) 

Since the field is degenerate and magnetic by (2.2), there exists a time- like frame at each 
point of space-time in which E vanishes (Komissarov 2001). This is known as a proper frame 
(not necessarily a coordinate frame). Denote the magnetic field in the proper frame as b. 

It is insightful to decompose a cold MHD plasma into two fiuids associated with the 
two species of charge and take the zero mass density limit to attain fide. The 4-current 
density is expressed in terms of the 4- velocities of the two fluids, = —en^u'^ + en^u't- 
The momentum equations of the two fluids are: 

d en± u+ X B\ , , 

Adding the two momentum equations and taking the limit of nm — > 0, yields the force- free 
condition (2.1). Notice that in the same limit, (2.3) implies that all particles flow parallel 
to b, otherwise the resulting E in the frame of the particles would drive large cross-fleld 
(nonforce-free) currents. 

The maximum velocity that a particle can attain in the force-free limit is the speed of 
light along b. Transforming this velocity from the proper frame to a general orthonormal 
frame by a local Lorentz boost E x B/i?^ yields the maximum three velocity of a force- free 
particle in a local frame Lightman et al (1975); Blandford (2002), 

V _ ^ExB + (B^-£5^ 

» max — ^2 ' V / 



3. The Charge Horizon 



It is useful to deflne the proper frame known as the corotating frame of a ffde magne- 
tosphere. For such a frame to be meaningful, the poloidal magnetic flux should be slowly 
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changing (so as to make the azimuthal electric field, negligibly small by Faraday's Law) 
as in the fiat spacetime, rotating, stationary magnetosphere of Levinson (2004). The corotat- 
ing frame is attained by a local azimuthal boost, v^, from an any orthonormal frame at fixed 
poloidal coordinate in an axisymmetric magnetosphere, Vp = /B^, where the superscript 
P represents the poloidal component. If there is a non- vanishing E'^ then a poloidal boost is 
required in addition to an azimuthal boost, which greatly complicates the algebra. At a light 
cylinder, E^ /B^ — ±1 by definition. The light cylinder is also an Alfven critical surface in 
ffde which is a one-way membrane for Alfven waves, no Alfven waves can cross the surface 
anti-directed to the bulk flow (Blandford 2002). Using (2.4), the poloidal component of the 
maximum velocity of a force-free particle is 



where the plus (minus) sign determines the direction of the bulk flow (i.e., it is positive in 
a pulsar wind and negative near the event horizon of a black hole). Outside of the light 
cylinder, | E^ / B^ \> 1 and inside the light cylinder | E^ / B^ |< 1- Consequently by (3.1), 
in the rotating magnetosphere of Levinson (2004) or in a pulsar, V^^^. = at the hght 
cyhnder, V^^^^ > inside the hght cylinder and V^^^. < outside of the light cyhnder. Thus, 
no charges or particles can cross the Alfven critical surface (light cylinder) anti-directed 
to the flow - the Alfven surface is a charge and particle horizon in an axisymmetric ffde 
magnetosphere. By (2.4), at the light cylinder, the anti-directed particles spiral azimuthally 
at the speed of light with no poloidal velocity. 



4. Black Hole Chcirge Horizons 

As expected by the equivalence principle, it is straightforward to show that the concept 
of a charge horizon can be extended to curved space-time by means of an explicit example. 
We consider the case of the magnetosphere of a rotating black hole described by the Kerr 
metric. There are two light cylinders and hence two Alfven critical surfaces in a black hole 
magnetosphere, one is the standard outer hght cylinder of pulsar physics that was discussed 
in section 3 and the other is an inner light cylinder that is associated with the dragging of 
inertial frames in the ergosphere (Blandford 2002). Consider the initial state of the simulation 
of Komissarov (2001). In ingoing Kerr-Schild (K-S, hereafter) coordinates, the initial state 
has only one component of F^^ that represents a purely radial magnetic field. Prom (2.3), 
as nm — >■ 0, one has F^^u^^^ = 0. Evaluating this in K-S coordinates, tells us that particles 
moving along b satisfy the condition d(f) = 0. A constant K-S azimuthal coordinate, d(f) = 
transforms to the constraint d4> = —{a/A)dr in Boyer-Lindquist (B-L, hereafter) coordinates. 
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where A = — 2 Mr + a^. Note that this imphes a strong toroidal magnetic field near the 
event horizon in B-L coordinates (Punsly and Bini 2004). The d0 = condition defines 
the particle trajectories that are restricted to the field lines, a requirement of ffde that was 
derived in section 2. As discussed in Punsly and Bini (2004), the inner light cylinder is 
at the outer boundary of the ergosphere, the stationary limit in the initial state (the initial 
state has zero field line angular velocity as viewed from asymptotic infinity, Qp =0). In 
order to evaluate (2.4) and (3.1), we need to introduce an orthonormal frame so that the 
concept of three- velocity is well-posed. The most famous such frame in the Kerr space-time 
is the ZAMO frame (see Lightman et al (1975), Punsly (2001) and references therein for 
a review of ZAMOs and the stationary limit surface). Using the ZAMO field values in the 
initial state, found in Punsly and Bini (2004), evaluated at the stationary limit (inner light 
cylinder) shows that V^^^^ = at the light cylinder, V^^^ < inside the light cylinder and 
^max ^ outside of the light cylinder as expected from section 3. 

One can get more insight into the Alfvcn critical condition by looking at the azimuthal 
velocity of a particle as measured in the ZAMO frames, c(3'^, that fiows outward (as viewed 
in the global B-L coordinates) inside of the inner Alfven surface. First transform the four 
velocity from B-L coordinates to the ZAMO frames, then apply this result to trajectories 
restricted to b, d(j) = 0, 

/ — -I \ 
: (4-la) 

(4-lb) 

■ O. dr 

where r is the proper time of the particle, the angular velocity of the ZAMOs as viewed 
from asymptotic infinity, Q = —g,j,t/g4>(p, is defined in terms of the metric in B-L coordinates 
and a = a/ A sin^ 0/ g4,4> is the lapse function (that describes the redshift of the ZAMO 
frames as viewed from asymptotic infinity) that vanishes at the event horizon. By (4.1), for 
outgoing trajectories, the maximum permissible azimuthal three-velocity (corresponding 
to dr/dt — 0) is equal to -c at the stationary hmit and decreases without bound as one 
approaches the event horizon. As in flat spacetime, an outgoing particle can spiral endlessly 
at the Alfven critical surface, but never flow through it. 

Consider the ramifications of (4.1c) to a hypothetical charge that appears to fiow out- 
ward as seen globally {dr/dt > 0) that is inside of the inner Alfven critical surface. Since 
there is a strong toroidal magnetic field extant near the horizon in the initial state, the 
physical charges flow on spiral trajectories. Inside of the Alfen critical surface, the resulting 
spiral flow velocity of a globally outgoing force-free charge is necessarily much larger than 
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the speed of light in order that the poloidal projection of the flow velocity be sufficient to 
escape black hole gravity. Thus any out-flowing charge in this region is acausal. 



5. Longitudinally Polarized Fast Wave 

According to eqns. (8) and (11) of Levinson (2004), a fast wave can propagate at the 
speed of hght with The causality question is whether there is a light-hke wave that can 
propagate changes in E\\. It is useful to look at abrupt discontinuities in order to resolve 
this question for two reasons. First of all, the simulations of Komissarov (2001, 2004) evolve 
by means of a Riemann solver that propagates discontinuities that make changes in the field 
variables. Secondly, the step wave analysis of MHD waves is a method that was specifically 
developed to provide clarity to the causal evolution of waves (since they are launched by a 
well defined piston) in the classic work of Kantrowitz and Petschek (1966). In the WKB 
notation of Levinson (2004), the oscillatory nature wave functions is given by 

gi(VV'(x)-x-tjt) _ g-ifc(x)(D^t-X„) (5-1) 

where is the fast wave phase speed and is a local coordinate along k. For non-dispersive 
light-like waves, the step function essentially contains information on all of the oscillatory 
modes since it is given by the Fourier composition, 

where e is an arbitrarily small positive number in the usual sense. 

If the light-like conjecture is physical then it must result from an analysis of abrupt 
discontinuities in cold MHD in the limit of zero plasma mass density. We will evaluate the 
discontinuities in the rest frame of the wavefront and take the limit of zero mass. Technically 
there is no light-like frame of reference, nevertheless our results will be well-defined in the 
limiting process for which time- like wavefronts exist. MHD discontinuities are solved for by 
considering the continuity of the stress-energy tensor, T'^'^, across the wavefront. Let x be 
the local normal coordinate to the wavefront and the upstream magnetic field, B is in the 
x-y plane and z lies in the wavefront surface. All we need consider to get the desired results 
is mass conservation, the frozen-in condition, which are 

nu^u° = constant , E -v x B = , (5-3) 
the antisymmetric components of Maxwell's equations 

Fap.j + F-ya-P + Fp-y.a — , (5-4) 
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and the continuity of one component of the stress-energy tensor, T^^, 

nnu'^u' - ■^{E'^E' + B^B') = , (5-5) 

in which aU of the quantities are evaluated downstream unless there is a subscript "u" and 
fjL is the specific enthalpy (note that by (5.3), i?^ = 0). The continuity of Bj. used in (5.5) 
follows from the V • B = condition in (5.4). A tremendous simplification occurs in (5.4) at 
the step wavefront, since to lowest order, all the singular terms must cancel (surface terms 
like delta functions), thus the normal covariant derivative does not depend on the connection 
coefficients and one has continuity of E^ and E^ at the wavefront, or 

dEy _dE^ _^ 
dx dx ' 

across the wavefront. Inserting (5.3) and (5.6) into (5.4), one gets 

written in terms of the field in the plasma rest frame, b. For riu = n, is the Alfven or 
intermediate wave speed in MHD (Punsly 2001). Taking the limit of zero mass density, (5.7) 
has two solutions, 

y = -2£ = c cos^ , or = , (5-8) 
K 

where 9 is the angle between b and the wave normal in the plasma rest frame and v is the 
force-free value of the Alfven speed in a proper frame (Punsly 2003). Thus by (5.8), a force- 
free discontinuity travels at the Alfven speed (which is not the speed of hght in general) , or 
it carries no E\\. 



(5-6) 



n 



+ 



x\2 



= , 



(5-7) 



6. Discussion 



It was shown in the last section that force-free fast waves cannot propagate changes in 
£■11, in agreement with Komissarov (2002). A discontinuity in E\\ requires a surface charge 
density, a — £'||/47r, by Gauss' law. There is no precursor to a and it is the source of the 
wave. In fact, it can be shown that the motion of a along the magnetic field line in the 

proper frame (the surface current) and Maxwell's equations entirely determine the fields in 
the downstream state (Punsly 2003). The Alfven surface in ffde is a one-way membrane for 
the physical charges that comprise a and therefore for waves carrying E\\ . Furthermore, the 
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existence of the charge horizon and the inherent charge of Alfven discontinuities provides a 
fundamental physical description of the Alfven critical surface in ffde. 

How do we explain the results of Levinson (2004)? That calculation was as attempt 
to collect higher order terms in a WKB approximation. The danger in doing this is that 
one must be rigorous in making sure that all terms of higher order are retained in all equa- 
tions from the beginning. The inconsistent equation is the normal component of Ampere's 
law, F^^.^ = AnJ^/c. It was not justified to ignore the transverse derivative, Vfr of the 
tangential magnetic field, Btan- From eqns. (8) and (11) of Levinson (2004), dE^/dt ~ 
{w / c){5E2/kR) = SE2/R in the notation of Levinson (2004) in which R is a cylindrical 
radius. Therefore dE-^/dt ~ Vtr x ^Bian (the change in the poloidal electric field in a 
force-free wind also results in a change in the toroidal magnetic field, in fact the two field 
components are approximately equal in the asymptotic field zone beyond the Alfven critical 
surface (Ogura and Kojima 2003)). In particular the normal component of Ampere's Law 
changes from dE^/dt = AtiJ^ in Levinson (2004) to Vtr x ^-Btan = ^t^J^ /c in the full 2-D 
discontinuity calculation. The failure of Levinson (2004) is that as the wavelength increases, 
plane waves are more and more inaccurate representations of the plasma waves in a rela- 
tivistic magnetosphere. One must find more exact 2-D solutions, with variations transverse 
to the propagation vector to learn more than the standard WKB approximation. The 2-D 
solutions should reveal any dispersive effects of the fast waves that result from boundary 
conditions as is the case in a vacuum waveguide or a plasma-filled waveguide (Punsly 2001). 
The axisymmetric 2-D fast wave discontinuities in an ffde black hole magnetospheres were 
determined in Punsly and Bini (2004). The oscillatory 2-D fast waves in a black hole mag- 
netosphere are linear combinations of spin-weighted spheroidal harmonics convolved with a 
radial function that is a solution of an extremely complicated differential equation (Teukolsky 
1973; Punsly 2001). 
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